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Abstract 

Discretizations of the Laplacian operator on non-hypercubical lattices are 
discussed in a systematic approach. It is shown that order a 2 errors always 
exist for discretizations involving only nearest neighbors. Among all lattices 
with the same density of lattice sites, the hypercubical lattices always have 
errors smaller than other lattices with the same number of spacetime dimen- 
sions. On the other hand, the four dimensional checkerboard lattice (also 
known as the Celmaster lattice) is the only lattice which is isotropic at order 
a 2 . Explicit forms of the discretized Laplacian operators on root lattices are 
presented, and different ways of eliminating order a 2 errors are discussed. 
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I. INTRODUCTION 



Lattice simulation is an invaluable tool in our studies of nonperturbative aspects of 
quantum field theories. By the discretization of the spacetime continuum, the original system 
with infinitely many degrees of freedom is approximated by one with only a large but finite 
number of degrees of freedom, which can be studied numerically. In practice, the lattice has 
almost always been chosen as hypercubical since the original paper by Wilson |IJ, although 
some "uncanonical" choices has also been discussed. Four- dimensional random lattices were 
discussed by Christ et. al. in Ref. |^-§J. The body-centered hypercubical lattices in four 
dimensions had been first studied by Celmaster f||-[7| and subsequently in fsHTlll- Almost at 
the same time four dimensional simplicial lattices were studied in |12|-[16|. (An interesting 
comparison between the last two uncanonical lattices can be found in Ref. flT7| .) 

Let's recall the definition of a lattice. A set of vectors {Bi} (1 < i < n) in m-dimensional 
space R m is linearly independent if 

n 

Y / x i B i = iff Xi = VI (1.1) 

i=l 

With such a set of linearly independent basis vectors {B{\ one can define a lattice as below. 
An n-dimensional lattice A n is a subset of R m (n < m) such that 

n 

A n = {X = Y J Z l B l :z l eZ] (1.2) 

i=l 

The lattice always contains the origin and is symmetric under reflection X — > —X in R m . 
Notice that we have not required that n = m as some lattices (like the triangular lattice 
in two dimensions) can be more conveniently represented when embedded in some higher 
dimensional spaces. Also note that, under this definition, random lattices are not lattices. 
Still, such a definition encompasses a huge set of lattices, including extremely skewed ones 
whose unitcells are long narrow parallelograms or parallelohedra. However, such highly 
asymmetric lattices are rarely used for simulations. Instead, physicists always use lattices 
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which are "highly symmetric" (we will give a mathematically precise meaning to this phrase 
later). In particular, there is a class of lattices called the "root lattices", which contains 
most of the "highly symmetric" lattices, at least in lower dimensions n < 4. While there is 
no real justification in neglecting the lattices which do not belong to this class, root lattices 
are so well studied that they provide a natural starting point for our investigation. As we 
will see, there are exactly three root lattices in four dimensions, namely the hypercubical, 
the body-centered hypercubical and the simplicial lattices^, i.e., exactly the cases which 
have been studied. 

It is important to compare the relative merits of these different lattices in simulations. 
Several papers mentioned above attempted to answer this question by simulating the same 
physical system on different lattices and comparing the results. These comparisons are very 
directly relevant to the practical applications of lattice simulations, but theoretically they are 
"unclean" in the sense that they are a priori sensitive to the physical problems in question, 
the actual algorithms used, and other non-geometrical factors. This paper attempts to study 
the same problem in a complementary way. We will obtain a general prescription on how to 
construct the Laplacian operator, the simplest possible differential operator, on any lattice. 
While the Laplacian operator is not directly relevant for lattice QCD, it is used in lattice 
simulations of strongly coupled scalar 4 theories, and it is logical to study the Laplacian 
case thoroughly before tackling the more complicated cases of the Dirac and the plaquette 
operators. By studying the discretization errors of the Laplacian operators mathematically, 
one can parameterize these errors in ways which depend only on the geometrical properties 
of the lattices. The questions we will be addressing include: 

• Are there any lattices on which the discretizations of the Laplacian operator on the 
nearest neighbors are free of 0(a 2 ) errors? 

• What are the merits and demerits of our canonical choice — the simple hypercubical 



1 This fact has been observed in Ref. [12|. 
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lattice? 

• What is the role of symmetry? 

• In general, how do we eliminate the 0(a 2 ) errors? (Answer: by involving sites which 
are not nearest neighbors.) Is it always possible? 

• In a given dimension, does the lattice which has the smallest 0(a 2 ) errors still work 
best after the 0(a 2 ) errors are cancelled? 

In this paper, we will try to provide answers to all these questions. The approach is 
"clean" in the sense that all the effects of algorithms or choices of physics problems are 
disentangled. On the other hand, it must be admitted that its relevance to practical lattice 
simulations is not immediate. Roughly speaking, this paper is more about lattice field theory 
than about lattice simulations, and the theorems obtained should be regarded as general 
"rules of thumb" instead of infallible laws in practice. 

This paper is organized as follows: In Section II we present a prescription (one shell 
discretization) for approximating the Laplacian on a lattice. Root lattices are defined in 
Section III, and the explicit forms of the discretized Laplacian operators on root lattices are 
presented in Section IV. Order a 2 errors are studied in Sections V and VI, while the issue 
of anisotropy is discussed in Section VII. In Section VIII a general prescription (two shell 
discretization) for eliminating the order a 2 errors is presented, and the order a 4 errors in 
this scheme are estimated in Section IX, followed by the conclusion in Section X. 

II. LAPLACIAN ON A LATTICE 

On hypercubical lattices, the discretization of the Laplacian operator is usually con- 
structed through the trapezoidal rule applied on each lattice direction. For example, the 
second derivative in the x direction of an arbitrary function at the origin is 

a 2 d 2 x (j)(0) = 4>{ax) - 20(0) + <j>(-ax) + ..., (2.1) 

and by summing up over all lattice directions, 
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a 2 d 2 (f)(0) = E (0(a£) - 20(°) + <K-ap)) +■■■ 

A 

= { E ( exp(a^) - 2 + exp(-a«9 M )) Wo) + . . . , (2.2) 

A 

where the /t's are the unit vectors along the lattice directions. 

The construction above depends so heavily on the existence of a set of mutually orthog- 
onal lattice vectors {/t} that it cannot be naturally generalized to non-hypercubical lattices. 
(Try the two dimensional triangular lattice, for instance.) To obtain a construction which 
is equally applicable to hypercubical and non-hypercubical lattices, we will make use of the 
coordinate-independent definition of the Laplacian operator. Recall that when satisfies 
the n-dimensional Laplace equation d 2 (f) = 0, the value of at any point is the average of 
its value on an (n — l)-dimensional hypersphere of arbitrary radius centered at the point in 
question. 

0(0) = / d// n _! <j>(X), (2.3) 
where is the unit measure on 5 n_1 . 

/ da n ^ = ( dQ^Q^-S = 1. (2.4) 

Jgn-l P " Jgn-1 " 2W 2 V ' 

So one can take the Laplacian d 2 cf) as an indicator of how much the equality Eq. (|2.3|) is 
violated. Just from dimensional analysis one knows that 

«9 2 0(O) = NljmQ d// n _i0(X)-0(O)Vr 2 , (2.5) 

where r is the radius of the hypersphere S 1 ™ -1 , and N is a proportionality constant one can 
easily determine to be 2n.| 2 j 

2 The value of ./V can be determined by considering the following electrostatic system satisfying 
the Poisson equation d 2 (j) = p in n dimensions. Consider a ball of constant charge density p and 
radius r. One can easily see that the potential difference between the surface and the center of the 
ball is pr 2 /2n. Substitute this back into Eq. ( p.5[) and ./V is found to be In. 
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On a lattice, however, is only defined on discrete lattice sites, and one can at best 
approximate the integral with the sum of over discrete points. It is useful to introduce 
some formalism here and define a shell S r of lattice sites with radius r > as the set of all 
sites with ||X|| = r. Obviously these <S r 's are non-empty only for discrete values of the r's, 
and the smallest of these values, p, is the separation between each lattice site and its nearest 
neighbors and is usually called the packing diameter of the lattice. We will call S p the first 
shell, and its cardinality, i.e., the number of nearest neighbors, denoted by r, is called the 
kissing number^. Now we approximate the integral of over S n by its average £ over the 
first shell 

£0(0) = ± <l>( X k), (2-6) 

where £ is defined as an operator over 0, which can also be expressed in terms of the partial 
derivatives. 

£<K0) = {i E exp(X fc -«9)}0(O). (2.7) 

x k es p 

The corresponding approximation of the Laplacian operator is 

D 2 0(O) = 2n(± £ 0(X fe )-0(O)y P 2 = 2n(£-l)0(O)/p 2 , (2.8) 

Xfr£S p 

The nonlocal operator D 2 defined above, called the one shell discretization of the Laplacian 
operator, is the main object of this study. How good is D 2 as an approximation of <9 2 ? Note 
that 

p^V(0)/2n = i £ 0(X fc )-0(O) = i £ (0(X fc )-0(O)) 

X k £S p X k £S p 

= E (exp(X fc -d)-l)}0(O). (2.9) 

X k £S p 



3 Also known to chemists as the contact number or the coordination number, and to historians as 
the Newton number, after Isaac Newton who studied the problem in three dimensions. 
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So the expression inside the braces is a derivative representation of the nonlocal operator 
p 2 D 2 /2n. By expanding the exponential, which is equivalent to a Taylor expansion in p, one 
can express p 2 D 2 /2n in terms of powers of the derivative operator. Now 

ex P (x fc • d) = i + x k ■ d + \{x k ■ d) 2 + \{x k ■ df + CV). (2.10) 

The order p° term is just unity, which is cancelled by the "—1" coming from the 0(0) term 
in Eq. (|2.8| ). The order p 1 term is proportional to 

( ]T X k )-d = 0, (2.11) 

x k es p 

due to the reflection symmetry of the lattice. Similarly, all odd powers of p vanish due to 
reflection symmetry. On the other hand, the p 2 term is | ^Z k {X k ■ d) 2 which is in general not 
a Laplacian unless the set of nearest neighbors S p has vanishing quadrapole moment. So we 
have the theorem: 

The one shell discretization D 2 of the Laplacian operator over the hrst shell S p is equal 
to the Laplacian operator d 2 itself in leading order if and only if the quadrapole moment of 
S p vanishes. 

It must be emphasized that the theorem above applies only for one shell discretizations as 
defined in Eq. (|2.8|) . In general one may give different weights to different lattice sites in S p 
in the sum in Eq. ( |2.8|) or even extend the sum to sites which are not nearest neighbors, and 
the theorem does not apply to such constructions. However, such alternative discretization 
schemes contain many free arbitrary parameters (as in the weights) and they are clearly 
beyond the scope of this paper. 

III. ROOT LATTICES 

We have seen that the one shell discretization method provides a good approximation of 
the Laplacian operator as long as the nearest shell of neighbors has vanishing quadrapole 
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moment. Unfortunately, the author is not aware of any classification of the lattices satisfying 
this condition. However, it is well known that root lattices have been completely classified, 
and all of them have vanishing quadrapole moment. Moreover, both the body-centered 
hypercubical and the simplicial lattices fall into this class, and so does the usual hypercubical 
lattice. So it is useful to study these root lattices closely. 

A root lattice is one whose basis vectors can all he chosen to be nearest neighbors and 

the angles between any two of them are either n/2 or 2n/3. In other words, when the basis 

vectors are scaled to have unit length, the inner products of any two of them are either or 
_i 

2 ' 

This definition depends solely on the relative directions of the nearest neighbors (as opposed 
to definitions in textbooks on Lie groups, which may involve next-to-nearest neighbors) is 
more useful for the study of one shell discretizations, which also just depend on the first 
shell. 

The classification theorem for root lattices is as follows: 

An n-dimensional root lattice is either isometric to one of the Cartan-Dynkin root lattices 
Z n , A n , D n and E n for some value of a, or is an outer product of several of them with the 
same a. 

The Cartan-Dynkin lattices Z n , A n and D n are defined as below for any n e Z: 
Z n = {X = {x 1 ,...,x n ) EaZ n }, 

A n = {X = (x 1 ,..., x n , x n+1 ) E aZ n+1 ; ]T x fc = °}' (3- 1 ) 

k 

D n = {X = (xi,...,x n ) G aZ n ; J2 x k G 2aZ}, 

k 

These representations have the nice features that the coordinates of the lattice sites are 
always integer multiples of the lattice spacing a, and each lattice is symmetric under per- 
mutations of the coordinates. It is obvious that Z n are the square, cubic and hypercubical 
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lattices for n = 2, 3 and 4 respectively. The simplicial lattices A n are for n = 2 the hexag- 
onal lattices, for n = 3 the face-centered cubic lattices, and for n = 4 the four dimensional 



simplicial lattices studied in Ref. |12|-|16||. The checkerboard lattices D n are also square lat- 
tices for n = 2 (D 2 — Z 2 ), face-centered cubic lattices for n = 3 (D 3 ~ v4 3 ) and Celmaster 
lattices 0-0] for n = 4. Lastly, the i? n lattices are defined only for n = 6, 7 or 8. Since 
we are mainly interested in lattices in low dimensions, we will not discuss E n and refer the 
interested reader to Ref. |1| instead. 

Several questions immediately come to mind. Firstly, are the A n and D n really the same 
Cartan-Dynkin root lattices one uses for the representations of Lie groupsQ? What is Z n in 
that language, and where are the B n , C n , and the exceptional groups? The first question 
can be easily answered by evaluating the angles between the basis vectors of A n _i. 

(x h ...,x n ) = (1,-1,0,..., 0),(0, 1,-1,0,..., 0),...,(0,..., 0,1,-1), (3.2) 

where each vector is rotated by an angle 27r/3 from its predecessor as described by the 
A n Dynkin diagrams. The D n lattice has one more basis vector than A n _i, namely 
(— 1, — 1, 0, . . . , 0), which is orthogonal to the first basis vector of A n _i but makes an angle of 
27r/3 with the second, exactly as decreed by the D n Dynkin diagrams. Lastly, the Z n lattice 
also has one more basis vector than A n _x, namely (— 1, 0, . . . , 0), which makes an angle 37r/4 
with the first basis vector of A n _i, is orthogonal to the rest, and is shorter than the rest of 
the basis vectors by a factor of 1/V2, i.e., it has exactly the properties of the "double bond" 
of the Dynkin diagram of B n . In other words, Z n is the root lattice of the Lie group B n . 

The last exercise illustrates the fact that a given root lattice can be spanned by different 
sets of basis vectors (and hence different Lie groups), and the choice most convenient for the 
study of Lie groups is not necessarily the most convenient for our purposes. As we have seen, 
the set of nearest neighbors S p plays a special role in our discussion. As a result, we have 



4 Readers who are not familiar with the Cartan-Dynkin classification scheme for Lie groups can 
consult Slansky and Georgi [ 20 1 . 

9 



chosen to define root lattices in a way which depends only on the relative directions of the 
nearest neighbors. The unique Lie group defined by the nearest neighbors on a given root 
lattice is always simply laced (no multiple bonds, all basis vectors have the same length). 
For A n it is SU(n + 1), for D n it is SO(2n) and for Z n it is SU(2) n (not SO(2n - 1)).Q 

The remaining Dynkin diagrams (C n , _F 4 and G?) are all multiply laced (with multiple 
bonds, and basis vectors of different lengths). One can easily check that the C n and D n 
groups have the same root lattices, but the former is doubly laced while the latter is simply 
laced. Similarly, the simply laced counterparts of G2 and F4 are A2 and -D4 respectively^. 

So we have clarified the relationship between the Lie groups (which have little to do with 
this paper) and the root lattices (our real concern). The reader should be warned that, 
for the purposes of this paper, it is advisable to temporarily forget about the Lie group 
structures of the lattices and regard the root lattices as geometrical constructions in their 
own right. 

IV. DISCRETIZED LAPLACIAN OPERATORS ON ROOT LATTICES 

In this section, we will give explicit formulae for E and D 2 on root lattices Z n , A n and 
D n . Recall that, since all lattices are symmetrical under reflection with respect to the origin 
(X m — > —X m for all m, m being the dimensionality of the coordinate space, i.e., n for Z n 
and Dn and n + 1 for A n ), in the Taylor expansion of E, only terms with even powers of p 
survive. Hence 

S-i^ exp(X fc • d) = O + a 2 2 + a 4 A + a 6 O e + . . . . (4.1) 

5 This is why we avoid calling the hypercubical lattice B n — to emphasize that we are studying 
the root lattice, not the Lie group. 
6 The authors of Ref. ipjdTll actually referred to the Celmaster lattice D4 as the F4 lattice. 
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The operators Oi are Z-th order derivatives, their forms constrained by the symmetries of the 
lattices. As mentioned before, all three root lattices defined in Eq. ( |3.1| ) are invariant under 
permutations of the coordinates. Moreover, Z n and D n are also symmetric under reflection 
of individual coordinates (X m — > — X m for a given m). These symmetry constraints limit 
the number of possible operators appearing at each order. 

To enumerate the number of possible operators appearing at each order, the following 
definitions are useful. 

^r'L, E <>!:<>!:■■■<>!:■ all fs even, 

E dildit some fs odd. (4.2) 

l<i 1 ^i 2 ^...^i p <m 

The definitions of the even operators B and the odd operators A are actually identical except 
that A may contain odd powers of <9;, which is not invariant under reflection of coordinate 
Xi. As a result, 0\ for Z n and £> n contains only terms like <Sj"j 2 ...j with X) P Jp — ^ but no A 
terms. On the other hand, for A n , both B^t}\ and A^t}\ contribute. 

J1J2--JP JlJ2---.Jp 

It is easy to see that Oq = 1 for any lattice. There are two operators which appear at 
order a 2 , namely 

— J2i<i^j< m didj. (4.3) 

As mentioned before, for Z n and D n , A does not contribute and O2 must be proportional to 
the even operator , which is exactly the Laplacian operator in m dimensions, and m = n 
for Z n and D n . On the other hand, for A n , O2 cannot be simply proportional to B^, as 
we want something proportional to the Laplacian operator in n, not m = n + 1, dimensions. 
However, for A n the odd operator A^ does contribute, and as we will see below, exactly 
in the right amount to give the Laplacian operator in an n dimensional subspace. 
When one goes to order a 4 , more operators appear. 

.o(m) _ „ 04 
°4 — 2^1<i<m °i 1 
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w(m) _ v -^2o2 
°22 — 2^1<i^j<m u i a j > 

Am) __ v a3a 
•"4-211 = J2l<i=/=j=/=k<m ^djdk, 

4m — T,i<i^j^k^i< m 9idjd k di. (4.4) 

As before, the odd operators only contribute in the v4 n case while the even ones in general 
may contribute for all three classes of root lattices. 

After discussing the generalities, we will move on to specific cases and calculate E, and 
hence D 2 = 2n(S — l)/p 2 , for all three classes of root lattices. We will also list the kissing 
numbers r, packing diameters p and densities of lattice sites 5 of the lattices. Many of the 



properties listed below are extracted from the excellent book by Conway and Sloane |L8 
although the definitions of p and 5 here are different from those in the book. 



A. The Hypercubical Lattice Z n 

The hypercubical lattice Z n , defined as 

Z n = {X = ( Xl ,...,x n )eaZ n }, (4.5) 

has kissing number r = 2n, packing diameter p = a and density 5 = l/a n . One can easily 
evaluate S, which takes the form 

£ = 1 + + i-yB^ + 0(a 6 ). (4.6) 

One can easy show that, due to the orthogonal nature of Z n , 

Oi = £ E # ( 4 - 7 ) 

l<i<n 

Lastly, the one shell discretization D 2 is 



D 2 = Bf + Ja 2 Bf + • • • = d 2 + ^a 2 £ df + 0(a 4 ), (4.8) 

l<i<n 

and as predicted, D 2 = d 2 in the leading order. 
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B. The Checkerboard Lattice D 



The checkerboard lattice D n , defined as 

D n = {X = ( Xl , ...,x n )E aZ n ; ]T x k € 2aZ}, (4.9) 

k 

has kissing number r = 2n(n— 1) and packing diameter p = \[2a. Only half of the lattice sites 
of Z n is occupied in D n , so its density 5 = l/2a n is half of that of Z n . It is straightforward 
to calculate S, 

S = 1 + i-yB { ^ + JLcMW + a^hya'fiW + 0(a 6 ), (4.10) 

and bearing in mind that now p 2 = 2a 2 , 

D 2 = Bt ] + |a 2 Bf + ^a 2 B$ + 0(a 4 ) 

= d 2 + ^- 2 a 2 df + j^a 2 £ d 2 d 2 + 0{a±). (4.11) 

Again, we verify that .D 2 = d 2 in the leading ordeiQ. 

The four dimensional checkerboard lattice is of special importance as it has been 
extensively studied by Celmaster JS|-0] and is usually known as the Celmaster, or body- 
centered hypercubical lattice. 

Dl = {X = ( Xl ,x 2 , x 3 , x 4 ) G aZ 4 U a(Z 4 + (|, ±, §, §))}, (4.12) 

i.e., the coordinates are either all integers or all half integers. The Celmaster lattice D\ 
is actually the dual, or reciprocal lattice, of D 4 . Reciprocal lattices are in general quite 
different: in three dimensions, the face-centered cubic and body-centered cubic lattices are 
reciprocal, and are not similar at all. In this special case, however, and D\ are actually 
similar to each other. First note that both have 24 sites in the first shell: 

D 4 : (±1, ±1,0,0) + perm.; 

D* 4 : (±1,0,0,0) + perm., (±±, ±\, ±\, ±\) + perm.; (4.13) 

7 This discretization is equivalent to that used in Ref. [pj-|Tl| 
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and 24 sites on the second shell 

£> 4 : (±2,0,0,0) + perm., (±1, ±1, ±1, ±1) + perm.; (4.14) 
Dl : (±1, ±1,0,0)+ perm. . 

It is clear that the first shell of D\ is exactly the second shell of D\ and the second shell of 
_D 4 is the first shell of D\ scaled up by a factor of 2. One can actually show that one can start 
out with D4, perform a global rotation and overall rescaling by a factor of l/\/2 and end 
up with Dl. So -D4 and Dl are geometrically similar in the sense that similar triangles are 
similar — they can be mapped onto each other through global rotation and overall scaling. 
As we will see below, the first and second shell of D4 being similar will be important when 
we try to cancel 0(a 2 ) errors by going beyond one shell discretization. 

C. The Simplicial Lattice A n 

Lastly, we will discuss the simplicial lattice A n , which is the hardest to study as it is 
embedded in an n + 1-dimensional coordinate space. 

A n = {X = (x u ...,x n ,x n+1 ) 6aZ n+1 ;^i fc = 0}. (4.15) 

k 

it has kissing number r = n(n+l), packing diameter p = \^2a, and density^ 5 = 1/ \/n + la n , 
and 

S = 1 + 2!(n+l) a2 ^2 + ) + 2!n(n+l) a2 ^ll + ) 

+ lK^I)^ n+1) + ^TI)^ +1) + IR^ 4 4i +1) + 0(a% (4.16) 

without any terms proportional to .4.211 or Ami- One can then read off D 2 , which takes a 
rather complicated form. 



8 The easiest way to show this is by observing that Z n +i has density l/a n+1 , and the ^4 n -like 
sublattices {X £ aZ n+1 ;J2k x k = C} where C is a fixed integer, have density \/n +T/a in Z ri+ \. 
The quotient of the two gives the density of A n . 
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_ (n+l) D 2 (n+D^H 



(n+1)' 

2n „2 »("+!) i 6n „2r}( n + 1 ) 6n 9 



(n+T) E ~~ (n+T) X! 

l<i<n+l l<«^i<n+l 



(4.17) 



+ T2#K)« 2 E ^ + i(in)« 2 E d 2 d 2 -^a 2 ^<^< n ^d) + o{a^. 

l<i<n+l l<«^j<n+l 

Does the sum of the first two terms give the Laplacian operator in the n-dimensional subspace 
described by J2 X % — 0? Note the following relation between the Laplacian operators in n 
and n+1 dimensions. 



R n + 1 



+ d 2 

~ u xi 



(4.18) 



where x — (1, 1, ... , l)/yn + 1 is the unit normal vector of the subspace J2%i — 0- O ne can 
easily check that 



dl=[ E 3* + 

l<i<n+l 



E 

l<i^j<n+l 



did^/(n + l), 



(4.19) 



and hence 



<9 2 



= 9 2 



R n+1 



dl = 



In \ I.) E ^ 

l<i<n+l 



(n+1) 



E 



l<i^j<n+\ 



(4.20) 



exactly as advertised. 



V. DISCRETIZATION ERRORS 

We have seen that D 2 does reproduce d 2 exactly in leading order for all three classes 
of root lattices we have studied. In this section, we will analyze the order a 2 errors of D 2 , 
which we will denote by E and which is related to O4 by 

E = 2n0 4 a 4 /p 2 . (5.1) 

Its forms for all three classes of root lattices are summarized below. 
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Z n : E — ^a 2 2B^\ 
D n :E= Ja 2 (2Bf° + 
A n :E = ia^^T 




+ 



n(n+l)"22 



6n »(n+l) 



n(n+l) 



6n 




(5.2) 



At first sight the coefficients may look random and without any discernible patterns, but as 
we will see below they are actually not as random as they look. 

All the expressions for E above are quartic polynomials of <9j. Naively, one may estimate 
their size by p 4 , where p is the typical momentum scale in question. However, such an 
estimate is sloppy, since as an operator p A = <9 4 = (<9 2 ) 2 , i.e., the square of the Laplacian 
operator, which is a scalar operator in n-dimensional Euclidean space, while B4, B22 and 
A31 are all non-scalar operators as defined. Thus motivated, we will define E , the scalar or 
isotropic part of E, as 



where d/i n -i is again the unit measure onS n 1 . In other words, E is E integrated over all 
directions in n-dimensional momentum space. When one evaluates E , it turns out that 



which is surprisingly simple (when contrasted with the expressions before the integration). 
Moreover, for all cases one has 



The above relation is not obvious (at least to the author) from the raw expressions for 
E, but once motivated it can actually be proved very easily for all lattices (not necessarily 
root lattices). Recall that £ is defined as 




(5.3) 




(5.4) 



E = p 2 «9 4 /4(n + 2). 



(5.5) 




(5.6) 



For each site X k e S p , one can expand the exponential as 
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exp(X fc • d) = 1 + ±(X k ■ d) 2 + ±(X k ■ df + odd terms + 0(p 6 ), (5.7) 

where the odd terms will eventually be cancelled and do not concern us here. The quadratic 
and quartic terms can be rearranged in the following way. 

(X k ■ d) 2 = p 2 d 2 /n + [(X k ■ d) 2 - p 2 d 2 /n 

(X k ■ <9) 4 = 3p 4 d 4 /n(n + 2) + [{X k ■ d) 4 - 3p 4 <9 4 /n(n + 2)j . (5.8) 
The integrals 

d Mn _ 1 (X k ■ d) 2 = p 2 d 2 /n, [ dp n ^ (X k ■ df = 3p 2 d 4 /n(n + 2), (5.9) 



decree that the terms in square brackets in Eqs. ( |5.8| ) vanish under integration over all 
directions. Hence 

/ cK-i E = \ E (1 + ^h n P 2 d 2 + lRt+2)/^ 4 ) + 0( P % (5.10) 

x k es p 

which is a sum over r identical expressions and then divided by r, and is equal to the 
expression itself, 

/ dp n . x E = 1 + ^- n p 2 d 2 + ^yp 4 9 4 + 0(p 6 ). (5.11) 

Lastly, we have the scalar part of D 2 . 

J dun-x D 2 = J d/i n _! 2n(S - l)/p 2 = d 2 + ^2)P 2 <9 4 + G(p 4 ), (5.12) 

and Eq can be read off simply as the p 2 term. 

The simple result Eq. ( |5.5|) has important implications. Firstly, it gives us the following 

impossibility theorem: 

One shell discretization of the Laplacian operator always suffers from order a 2 errors, 
regardless of the choice of the lattice. 

One might have hope that, by using a lattice with a high degree of symmetry, all 0(a 2 ) 
terms may vanish identically. We now see that is an impossible dream. There is always a 
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scalar part which is not prohibited by any symmetry. The only way to eliminate order a 2 
errors is by involving lattice sites which are not nearest neighbors, as we will discuss later 
in this paper. 

Another corollary of Eq. ( |5.5| ) is that the scalar discretization error Eq depends only on 
a single geometrical parameter — the packing diameter p. (It also depends on the spacetime 
dimensionality n, but that is usually determined by the physics problem in question.) One 
can then ask, in a given spacetime dimension, which lattice will minimize Eq? It will be the 
subject of the next section. 



VI. ERROR DENSITIES OF ONE SHELL DISCRETIZATIONS 

As defined, the scalar discretization error E is independent of how we describe the lattice 
on the coordinate space. This is certainly a pleasing feature to physicists, who know that 
all physical quantities should not depend on the choice of coordinates. On the other hand, 
Eq is proportional to p 2 , which is not surprising; the smaller p 2 , the more lattice sites we 
are using to describe a given physical volume, and hence the smaller is the error. In other 
words, the error grows with 1/5, the volume of the unitcell of the lattice. To determine 
which class (Z n , D n or A n ) is optimal for our purposes, one defines the error density e as 
the ratio of the E to (a certain power of) the unitcell volume. 

The —2/n power has been chosen such that e is p independent; recall that E Q and 5 scale 
like p 2 and p~ n respectively. Hence e thus defined is a sheer parameter (not an operator) 
which just depends on which lattice which are using. For root lattices, it carries the value 

f 1, for Z n \ 

2 l-2/n for Dn . ( 6 2 ) 

t 2(n + l)~ 1/n , for A n . 
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Obviously one can check that the error density of Z n is always smaller than that of D n and 
A n .f\ This is certainly not a surprising result, as we have shown in the previous section that 
Eq is proportional just to p 2 , and for our root lattices with the same p, Z n is the one with 
the smallest density. 

We have seen that, among root lattices, Z n has the smallest error density. We can 
actually prove the much stronger result that Z n has the smallest error density among all 
lattices whose S p have vanishing quadrapole moment, i.e., all the lattices with D 2 = d 2 in 
leading order. The proof runs in two steps. First, observe that among all lattices whose 
basis vectors can be chosen to be nearest neighbors, Z n has the lowest density. (Among all 
parallelohedra with all sides of unit length, the one with the largest volume is a cube.) So Z n 
has the lowest error density among them. What about lattices whose basis vectors cannot 
all be chosen as nearest neighbors? Note that, for such a lattice A, there exists a proper 
sublattice Ac A which is generated by the nearest neighbors of A. The lattices A and A 
have the same p and hence the same Eq, yet being a subset, the density of A is smaller than 
that of A. As a result, the error density of A is larger than that of A, which in turn by the 
previous argument is larger than that of Z n . So we have the important result: 

Under one shell discretization, the hypercubical lattices have the smallest error densities 
among all possible choices of lattices in the same number of spacetime dimensions. 

Most lattice simulations are performed on hypercubical lattices, but usually the reasons 
are historical (following Wilson's example in Ref. |IJ) or algorithmic (D 2 can be easily written 
down as in Eq. ( |2.2| ) by using the trapezoidal rule). Here we provide another rationale: the 
hypercubical lattices have the lowest error density. 



9 Moreover, one also has ez 2 = £d 2 an d e^ 3 = €d 3 , i-e., the error density of a lattice is independent 
of how it is named. 
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VII. ANISOTROPY 



While Eq measures the isotropic part of the order a 2 errors, E — Eq measures the 
anisotropy of D 2 . The operator E — E is a quartic polynomial in <9j, which in general 
can have all £-th multipole moment up to I = 4. However, the first and third (dipole and 
octupole) multipole moments are identically zero as demanded by reflection symmetry, the 
second (quadrapole) moment is explicitly required to vanish, and the zeroth (monopole) 
moment Eq has been subtracted out. So E — E Q is purely of the fourth (hecadexapole) 
moment. 

It is straightforward to verify that E — E is non-vanishing for all root lattices with one 
exception: £) 4 , the checkerboard lattice in four dimensions. 

E = pa 2 (B 4 + B 22 ) = la 2 d 4 = E for D A . (7.1) 

That D 4 is more isotropic than Z 4 has been noted by Celmaster in Ref. H], but here the 
result is stronger. 

The checkerboard lattice in four dimensions D± is exactly isotropic at order a 2 . It is the 
only unexceptional root lattice with this property. 

One can understand the isotropy of the D4 lattice by studying its geometry. The D± lattice 
has an accidental threefold discrete symmetry (which is also a symmetry of its Dynkin 
diagram) which mixes £>4 and B 2 2- The only combination which is invariant under this 
threefold symmetry is <9 4 . In other words, rotational symmetry is protected by this accidental 
threefold discrete symmetry at order a 2 . 

Since we are mainly interested in lattices in low dimensions (n < 4), we have not studied 
the anisotropy of E 6 j )8 , although the calculation should be straightforward (though tedious 
- the packing number r is 240 for E 8 ). We also have not studied the issue of anisotropy 
for non-root lattices. However, since we have seen that the isotropy of -D4 is related to an 
accidental discrete symmetry, a feature not apparent in these other cases, the author would 
be very surprised if other lattices share this property of isotropy. 
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VIII. BEYOND ONE SHELL DISCRETIZATION 



We have seen that order a 2 errors are unavoidable in one shell discretization. The order 
a 2 errors can be substantial on coarse lattices; for example, as a rule of thumb the order 
a 2 error on a hypercubical lattice is around 10% if a is chosen to be around a third of the 
size of the physical system in question ||21[|. In order to reduce the error, one may try to 



use finer lattices. It is, however, a very costly procedure as the cost of lattice simulations 
grows at least as fast as the number of lattice sites, which is proportional to a~ n . (And 
that is before taking critical slowing down into account.) For n = 4, to reduce the error 
to 1% would increase the cost a hundredfold. Another way of seeing the same result is by 
observing that the error density e, which is a independent, is the proportionality constant 
between E and 5~ 2 ^ n . To reduce E by a factor of 10, 5 must be decreased by a factor of 
10™/ 2 , i.e., a hundredfold in four dimensions. 

Clearly it is not practical to control the a 2 errors just by using finer and finer lattices. 
Instead, one can completely eliminate the a 2 error by involving sites not in S p , i.e., by 
going beyond one shell discretization. As an example, we will review how this is done on 
hypercubical lattices. 

Recall that, for hypercubical lattices, one shell discretization on S a gives 



D 2 [S a ] = Bi n) + %a 2 Bt ] + O(a') = d 2 + ±a 2 £ df + 0{a% (8.1) 

Ki<n 



The argument in the square brackets is to remind us of which shell we are using. One can 
also discretize on the following set of 2n sites: 

2S a = {{±2a, 0, . . . , 0) + perm.}. (8.2) 

Notice that 2S a is just S a scaled up by a factor of 2. As a result, 

D 2 [2S a ] = 8 2 + ^(2a) 2 £ dt + 0{a A ). (8.3) 

\<i<n 

Consequently, the combination |-D 2 [iS a ] — ^D 2 [2S a ] is free of order a 2 errors. 
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Notice that 2S a is not the second shell of Z n . (In fact it is just part of the fourth shell.) 
Since the errors are in general proportional to some high powers of the radii of the shells 
involved, can we do better by using just the first and second shell, S a and «S V ^ ? Recall that 
^V2a i s the first shell of D n , and D 2 has been calculated to be 

D 2 [S^ a \ = Bi n) + |a 2 4 n) + ^a 2 B& + 0(a 4 ), (8.4) 

which contains both B^ and B22 . In contrast, D 2 [S a ] involves just the former but not the 
latter. Thus one cannot obtain an expression free of both B^f 1 and £> 22 ^ errors by taking a 
linear combination of D 2 [S a ] and D 2 [S^2 a ]. In other words, one cannot eliminate the order 
a 2 order by taking a linear combination of D 2 of any two shells. 

Bearing the above examples on the hypercubical lattices in mind, we will now move on to 
discuss the framework of elimination of the order a 2 errors on lattices in general. Consider 
the one shell discretization on two shells of radii r\ )2 with r\ < r 2 . Then we want to see if 
the expression 

D 2 [S ri ,S r2 ] = Cl D 2 [S ri ] - c 2 D 2 [S r2 ], (8.5) 

is equal to the Laplacian operator without any a 2 errors with some choices of coefficients 
ci,2- Obviously in order to get the coefficient of the Laplacian operator to be unity we need 
Ci — c 2 = 1. We also want the scalar errors E of the two shells, which are proportional to 
r 2 , to cancel, by demanding c\r\ — c 2 r\ = 0, or Ci/c 2 = (7"i/r 2 )~ 2 . Putting these together 
we have 

ci = r\l{r\ - r 2 ), c 2 = r 2 /(r 2 - r\). (8.6) 

So now we can define the two shell discretization of the Laplacian operator over shells 
S ri and S r2 as 

D 2 [S ri ,S r2 ] = (rlD 2 [S ri \-rlD 2 [S r2 \}/{rl-rl) 

= d 2 + a 2 H 2 + a' l n 4 + 0{a 6 ). (8.7) 
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Thus denned, the scalar discretization error of order a 2 is guaranteed to vanish, i.e., 
J d/j n -i TC2 = 0. However, 7i 2 itself may be non-zero, as shown by the example with 
the first and second shell in Z n . 

When we defined one shell discretization early in this paper, we had to explicitly demand 
that the shell in question have zero quadrapole moment. Here the situation is similar. For 
two shell discretization we have to demand that the two shell configuration in question, after 
weighting by ci t 2, has vanishing hexadecapole moment. In other words, the hexadecapole 
moments of the two shells must cancel each other. 

In principle, one can translate the condition of the cancellation of hexadecapole moments 
into mathematical constraints and look for their most general solutions — but in practice it 
is almost impossible. (We have not even done that for the quadrapole case!) For example, 
for Z n or D n , the shells described by (y±, . . . y n ) and (zi, . . . z n ) (plus permutation) will have 
their fourth multipole moment cancelled if and only if 

J2vI = KJ24, and E^ = ^E^ for some K. (8.8) 

i i i,j i,j 

These are coupled diophantine-like constraints, which in general cannot be solved. The 
corresponding constraints for A n are even more hideous. 

However, there are special cases which we know the hexadecapole moment are guaranteed 
to cancel, namely when the two shells are identical up to scaling (but no rotation). To be 
exact, let's define, for any natural number p > 1, 

pS p = {pX:Xe S p }. (8.9) 

Then the two shell discretization D 2 [S p ,pS p ] has vanishing hexadecapole moment and hence 
is free of order a 2 errors. The proof is immediate. For D 2 [S P ] = d 2 + p 2 Ti,2 + C*(a 4 ), one has 

D 2 [S p ,pS p ] = ({pp) 2 D 2 [S p ] - P 2 D 2 [pS p ]) I ((pp) 2 - p 2 ) 

= (p 2 P 2 (d 2 + p 2 n 2 ) - P \d 2 + p 2 p 2 n 2 )) /(p 2 - i) P 2 + o( P 4 ) 

= d 2 + 0(p 4 ). (8.10) 
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In practice, the two shells should be chosen as small as possible to minimize the error. 
This suggests setting p = 2, and 

D 2 [S P ,2S P ] = ±D 2 [S P ] - \D 2 [2S P \. (8.11) 

On the other hand, there is a unique case where the hexadecapole moments do cancel for 
two shells which are not related simply by scaling. As the reader may have guessed, this is 
the case for the first and second shells of the four dimensional checkerboard lattices, where 
the hexadecapole moment of each shell vanishes identically as dictated by the threefold 
discrete symmetry. Recall that both the first shell S^2~ a and the second shell S 2a contain 24 
sites as described in Eqs. ( [4.13 ,4.14), and the two shell discretization over these two shells 



is 

D 2 ^, S 2a ] = 2D 2 [S v - 2a ] - D 2 [S 2a ] for D 4l (8.12) 
which can be checked to be free of order a 2 errors. 

IX. ERROR DENSITIES OF TWO SHELL DISCRETIZATIONS 

Consider a two shell discretization with a 2 errors completely cancelled. 

D 2 [S ri ,S r2 ] = d 2 + a A H A + ... (9.1) 

To study the order a 4 errors, one needs to know 7^4, which in general is a sextic polynomial 
of <9j. To calculate TL± would be a painful exercise. However, it is easy to evaluate E' , the 
scalar part of TC4, 

E' = J dn n -i a A H A . (9.2) 

Again we will start with the Taylor expansion of exp(Xfc • d). 

exp(X fc • d) = 1 + ±{X k ■ d) 2 + ±{X k ■ 8f + ±{X k - 8f + odd terms + C(a 8 ). (9.3) 

With Xfr G S r (not necessarily the first shell), it can be rewritten as 
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exp(X fe • d) = 1 + \ { l 2 r 2 d 2 + \ x hr^ + ±I 6 r 6 d 6 

+ [...]+ odd terms + 0(r 8 ), (9.4) 

where l\ is defined by 

J dim-x{X h -d) 1 = hr l d\ (9.5) 

with I 2 = 1/n and I4 = 3/n(n + 2). We have not calculated Iq which fortunately is not 
needed for our purposes. Lastly, the [. . .] in Eq. ( |9.4D , like the terms in square brackets in 
Eq. ( [5.8|) , vanishes under integration over all directions. 
The scalar part of the one shell discretization on S r is 

J dfin-! D 2 [S r ] =d 2 + %( h /I 2 y 2 d* + !(/ 6 // 2 )r 4 <9 6 + 0(r 6 ). (9.6) 

and the scalar part of the two shell discretization on S ri and S r2 is 

D 2 [S r ,,S r ,] = d 2 - §(I e /I 2 )r 2 r 2 d G + 0(r 6 ), (9.7) 

and hence 

E>=-%(I 6 /I 2 )r 2 r 2 d«. (9.8) 

Now one can define the error density e' of a two shell discretization as 

^ = -|(/ 6 // 2 )e'5- 4 /^ 6 . (9.9) 

As expected, using larger shells will increase the error density. For root lattices, and using 
shells S p and 2S P , the error densities are 

' 4, for Z n ; 

e > = J 2 4 " 4 / n , for D n ; (9.10) 

t 16(n+ iy 2/n , for An. 
Again, the error densities are the smallest for hypercubical lattices. The reason is the same 
as the one shell case: for the same p, Z n has the lowest density. And the result can be as 
easily generalized to non-root lattices as before. 
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On the other hand, for D 4 we can simply use the first and second shells, and the error 
density is 

e = 4, for D 2 [S v - 2a ,S 2a ], (9.11) 
just as low as that of Z n . Hence we come to the following conjecture: 

Under two shell discretization, the hypercubical lattices have the smallest error densities 
among all possible choices of lattices in the same number of spacetime dimensions. The only 
other case with an error density as low is two shell discretization on the first and second 
shells of the four dimensional checkerboard lattice. 

The statement above is just a conjecture, not a theorem, as there may be two shells configu- 
rations other than S p and 2S P which have even lower error densities, just as the _D 4 case we 
have seen. However, exact cancellations of the hexadecapole moments of two shells are very 
unlikely to happen without a symmetry rationale (like the threefold discrete symmetry of 
D4), and since we do not see such extra symmetries in other lattices, the author concludes 
that it is very unlikely that counterexamples to the above conjecture will be found. 

X. CONCLUSION 

Let's recapitulate what we have learned through this study: 

• One shell discretization faithfully reproduces the Laplacian operator (with order a 2 
errors) if and only if the shell has vanishing quadrapole moment. 

• The order a 2 error for one shell discretization is non-zero for any choice of lattice. The 
hypercubical lattices have the smallest error densities in any spacetime dimensions. 

• Due to its threefold discrete symmetry, the four dimensional checkerboard lattice is 
isotropic up to order a 2 . It is the only root lattice, and also probably the only lattice (root 
or non-root), with this property. 

• Two shell discretization reproduces the Laplacian operator without order a 2 errors if 
and only if the hexadecapole moments of the two shells cancel each other. 
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• The order a 2 error for one shell discretization is non-zero for any choice of lattice. It is 
likely (though not yet proved) that the hypercubical lattices have the smallest error densities 
in any spacetime dimensions. In four dimensions, the error density of the checkerboard lattice 
is just as small as the hypercubical. 

One can see a clear pattern recurring throughout this discussion. The scalar (isotropic) 
properties of discretizations can be very cleanly studied by Taylor expanding exp(Xk-d) and 
then integrating over all directions. By involving more shells, one can eliminate the errors 
order by order. On the other hand, the anisotropic properties are in general very difficult 
to study. The isotropy of £> 4 is more like an accident (at least in our framework) — and 
it is straightforward to show that that the symmetry does not prevent higher multipoles 
appearing at order a 4 . As a result, the most tractable way of eliminating anisotropic errors 
is by cancellation between shells related just by a simple scaling. 

We have discussed the root lattices in detail and calculated their E explicitly. However, 
most of our results are equally valid for root and non-root lattices. However, root lattices 
are so tractable in theory, and useful in practice, that it makes a very natural starting 
point. The only important non-root lattice in low dimensions is D 3 , the three dimensional 
body-centered cubic lattice, which is the dual of D 3 . 

D* 3 = {X = ( Xl ,x 2 , x 3 ) G aZ 3 U a(Z 3 + (±, §, §))}, (10.1) 

Its coordinates are either all integers or all half integers. Its nearest neighbors are 
(±|, ±tj, ±|), and angles between them are not rational fractions of n. One can study 
.Dg just as we study the root lattices, but we do not find any properties of interest. 

What are the implications of this study for real lattice simulations? In no way does 
it challenge the canonical choice — the hypercubical lattice, which always has the smallest 
error density. The only case when one may consider otherwise is when isotropy is important; 
in such scenarios the four dimensional checkerboard lattice may be a tempting alternative. 
Still, as disclaimed in the introduction, the geometry of the lattice is only one of many factors 
which affect the efficiency of lattice simulations. Algorithmic considerations may overwhelm 
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geometrical ones; remember the two shell discretization of D 4 involves 48 neighbors, opposed 
to 16 for Z4. Still, the author expects this study is useful in providing a general framework 
in which different lattices may be discussed on equal footing, and the results may serve as 
"rules of thumb" for lattice simulations. 

In this paper, we have studied how the geometries of non-hypercubical lattices affect 
the discretization of the Laplacian operator. Possible generalization to the case of the Dirac 
operator is possible but not straightforward, due to the problem of parity doubling. Without 
doing a careful analysis, it is difficult to foresee whether the results we have obtained in this 
paper can be extended to the Dirac case, However, the author is willing to make the following 
conjecture: the order a 2 discretization error of the Dirac operator does not vanish for any 
choice of lattice. Just as the discretization of the Laplacian operator is always plagued 
with errors like a 2 <9 4 which have the same tranformation properties as the Laplacian d 2 , it 
is probable that the discretization of the Dirac operator 7 • d will suffer from errors like 
a 2 (7 ■ d)d 2 . Hopefully the author can extend the present framework to the Dirac case in the 
near future and determine the correctness of this conjecture. 
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